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Abstract 

66 ■ 

■ Recently, by analyzing the measurement data of Nikuradze, it has been proposed (N. Goldenfeld, 

o ■ 

CN ; Phys. Rev. Lett. 96, 044503, 2006) that the friction factor, /, of rough pipe flow obeys a scaling 

U ■ 

law in the turbulent regime. Here, we provide a phenomenological scaling argument to explain 



this law and demonstrate how intermittency modifies the scaling form, thereby relating / to the 
intermittency exponent, rj. By statistically analyzing the measurement data of /, we infer a 



Ch i satisfactory estimate for r] (~ 0.02), the inclusion of which is shown to improve the data-collapse 

' curve. This provides empirical evidence for intermittency other than the direct measurement of 

^ ■ velocity fluctuations. 
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I. INTRODUCTION 



A major aspect of fully developed turbulence is the existence of universal scaling laws in 
the so called inertial-range scales, ^ ^ ^ ^d, k {h) being the length scale of the energetic 
(dissipative)-range eddies. In particular, the scaling behavior of the 2nd-order velocity 
structure function under the Kolmogorov refined similarity hypothesis |l| is expressed as 

(1) 

where {si) is the average dissipation rate over a sphere of size / and r] is the intermittency 
exponent, whose value is considered universal. 

An important physical quantity in rough pipe flow is the friction factor /, which is related 
to the pressure drop across the pipe according to the Darcy-Weisbach formula (see e.g. j^). 
The friction factor is a function of the Reynolds number, i?e, of the flow and the relative 
roughness, of the pipe (r being the average size of the roughness elements and, R, the 
radius of the pipe). / can be expressed in terms of the wall stress, r, as / = where p 
is the density and V is the cross-sectional average of the mean-time-average velocity of the 
flow. 

In a seminal series of experiments on rough pipe turbulence that has remained a bench- 
mark in the fleld, Nikuradze [s] elucidated how f{^,Re) depends on its arguments. His 
data was presented as six curves, which are shown in Fig. 1. The main features of these 
plots are as follows. Up to Re ~ 3300, the flow is still laminar and / ~ This, of 
course, corroborates the exact result f = obtained from the Navier-Stokes equation. At 
Re ~ 3500, the entire curves plunge to the "smooth-pipe zone" in accord with Blasius's scal- 
ing 4, ^, / ~ Re~^ . Here the friction factor is independent of boundary roughness, due to 
the existence of a sufficiently thick viscous sublayer near the wall that screens the roughness 
elements from the turbulent flow. As Re increases, the roughness elements become progres- 
sively exposed to the turbulent flow and at large enough Re, the flow enters the "rough-pipe 
zone" , where the friction factor becomes dependent only on boundary roughness in accord 



with Strickler's law [6], / ~ (4; Is 



n 

R' 



These general characteristics can also be seen in other pipe [3, 1^, Is] and open-channel 
flow data. 

Recently, Goldenfeld has pointed out HI] that Nikuradze's data conform to the scaling 
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FIG. 1: Friction factor, /, of rough pipe flow versus Reynols number, Re, for different values of 
relative roughness, in Nikuradze's experiment. 



form 



f = Re Rei] 



(2) 



where the scahng function g{x) has the asymptotic behavior, 



9[X) 




and data collapse for fRe^ versus -^Re^ occurs for data that lie between the Blasius and 
Strickler regimes. The residual deviation of data collapse from the scaling form ([2]), apart 
from uncertainties in the data, may also reflect something more fundamental HI] . Here, we 
address the question as to what extent (if any) intermittency corrections, which have not 
been included in the original analysis, improve the data-collapse curve. 

Although turbulent flow in a pipe is anisotropic and inhomogeneous, Kolmogorov's theory, 
which rest's on the assumptions of isotropy and homogeneity, still applies l^, 14, Q, 
IGj. Here, using ([T]), we provide a scaling argument to explain the scaling form ([2]) and 
demonstrate how it is modified by intermittency, thereby relating / to r]. The measurement 
data of / can be, thence, utilized to infer the numerical value of rj, the inclusion of which 
is shown to improve the data-collapse curve. By statistically analyzing Nikuradze's data, 



we thus obtain an estimate for r] 0.02) that reconciles with the available values 



18 



19 



m 



2l\ . This provides empirical evidence for intermittency other than the direct 



measurement of velocity fluctuations. 



II. PHENOMENOLOGICAL SCALING ARGUMENT 

The shear stress exerted by the flow on the wall of the pipe scales as [12|, Il3| r ~ pVvr, 



where Vr = y {^Vr) is the velocity of eddies of length scale r, the size of roug hness elements. 
(Such eddies have dominant role in the momentum transfer to the wall |l2l. Il3l|). Thus, the 
friction factor can be written as / ~ y and hence from ([T]), / ~ where a = | + |. Let 
us now rescale such that r — Ir. Under this rescaling, we have / l°'f. Moreover, for 
Z > 1 (/ < 1), the rescaling results in a mere amplification (de-amplification) of roughness 
elements, without changing their geometry, thereby rendering the roughness to be exposed 
to the turbulent flow at lower (higher) Re. Thus, under the rescaling r Ir, we expect the 
same flow, albeit with Re l~^Re, where /3 > 0. That is, under rescaling. 

Collecting results, we therefore have 

f{lLj-PRe) = rf(L^Re) (3) 

Equation ([3]) indicates that / is a generalized homogeneous function of its arguments in 
the inertial range and the behavior is self-affine. We can expound further on the exponent 
(3, if we bear in mind that we are considering boundary layer turbulence. The turbulent 
transfer of momentum to the boundary wall by eddies, which predominantly takes place 
on length scale r (the size of dominant eddies), gives rise to an extra (eddy) viscosity near 
the wall. The effective viscosity in the boundary layer, therefore, scales as ~ rvr (length x 
velocity), i.e. as ~ r"+^. Hence on rescaling. Re l^'^^^Re, i.e. /3 = a + l = | + | and ([3]) 
reads: 

f(lLj-»-iRe) = l--f(L^Re) (4) 

Goldenfeld's scaling form ([2]) is just a particular form of (jl]), which can be obtained by taking 
/ = i?e"+T and ?7 = 0, of course. In presence of intermittency, the scaling form thus becomes 

2+3?) r 6 

f = Re s+3v g[— Re>i+3v) (5) 
R 



where, asymptotically, 

{const., X ^ 

Note that the Blasius and Strickler laws are, therefore, modified according to 



1 1 "7 

2+3i? / r \ 3"'"2 

f ^ Re 8+3. , / ^ i^-j 

respectively. The above modified Strickler formula coincides (to lowest order) with the result 
obtained by Gioia et. al. [l^. The modified Blasius formula, however, does not. The source 

3 

of discrepancy lies in their application of the expression Re~i for the Kolmogorov scale, 
which is invalid when the intermittency exponent is non-zero. The correct expression is 

s 

Re 8+3') , the application of which would have had resulted in the same modified Blasius 
formula as given above ^|. 

Scaling form ([5]) relates the friction factor to the intermittency exponent and provides 
ground for the estimation of the latter via empirical data other than the direct measurement 
of velocity fiuctuations. 



III. ESTIMATION OF THE INTERMITTENCY EXPONENT 

Having established (I5l), we proceed to infer the value of rj from the measurement data of 

2+3i7 6 

/. To this end, we plot fRe^+^'^ versus ^Re^+''^^ for different trial values of r] and examine 



the resulting data collapse by applying regression analysis (see e.g. |23l]). The best value for 
r] is one that will result in the highest possible data correlation, i.e., the smallest possible 
data scatter (deviations from data collapse curve). The data correlation coefficient, p, is 
defined by the fraction of the total variation (sum of squares) of data, y, that is explained 
by the regression model (fitting curve), according to 



ny-yf 



E{y-yy' 

Here y represents the average of data values, y, and y are their estimated values based on 
the regression curve. The closer p to 1, the more correlated the data, and the better the 
regression is. For a given regression model, we therefore examine the variation of p with 
rj] the best estimate of i] is the one that yields the correlation coefficient closest to 1. Now 
polynomial regression is suitable for nonlinear data curves that show a maximum/minimum 
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and the absence of this feature in the data collapse curve, as shown in Fig. 2 of reference [11 1 
(a non zero value of 1] does not change this feature), renders the polynomial model irrelevant. 
We, therefore, resort to the exponential regression model, the result of which is presented 
in Fig. 2. The estimated value of r] thus obtained is 0.02, which agrees with the gene rally 



■ees witn t ne g ene 

Q, 0, Q, HE 



expected range of values (~ 0.02 — 0.03) reported in the literature Il|,ll7|,ll8|,ll9|,l20|,l2l|. 
That is to say, Nikuradzes data set best conforms to r] ^ 0.02 and not, in particular, to 
r] = 0. It is rather interesting that large-scale properties like friction factors can provide 
evidence for intermittency, which is a direct manifestation of the small-scale statistics. This 



and similar previous observations [2^ embody deep connections between spectral structure 



and the global properties of turbulent systems. 
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FIG. 2: (Color online.) Exponential regression model. Best data correlation is obtained for the 
value 7] = 0.02 (top), resulting in the best data collapse (bottom). 
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